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The effects of resonances and flow on the correlation function for two identical particles 
are described assuming chaotic sources and classical propagation of particles. Expanding to 
second order in relative momenta, the source sizes can be calculated directly and understood as 
contributions from various fluctuations in the source. Specific calculations of source size radii 
are given assuming Bjorken longitudinal flow with additional transverse expansion. Results 
are compared to recent hit and KK correlation data from relativistic nuclear collisions with 
particular attention to the reduction in the correlation function due to resonances and the 
decreasing source sizes with increasing transverse momenta of the particles. 
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I. INTRODUCTION 

The HBT effect based on interference of identical particles is an important method in relativistic heavy 
ion collisions for extracting information about the source sizes and life times. || The correlations due to 
. interference are well described for pp collisions Q when the effect of resonances in particular the long lived 
' ones are included. In high energy nuclear collisions we hope to extract the spatial, temporal and momentum 
t-H , distribution of particles in the source at decoupling or freeze-out by correcting for resonances, Coulomb effects, 
and possible other final state interactions. The source sizes, life times and flow effects and their dependence 
on collision energy, impact parameter, projectile and target mass will be crucial for determining whether a 
quark-gluon plasma is created. Recent data allow determination of longitudinal, outwards and sidewards 
source radii for pions and kaons in heavy ion collisions at CERN and Brookhaven energies. The rapidity and 
transverse momentum dependence of source sizes have also been measured. ]5],^| The rapidity dependence of 
the longitudinal source size agrees with longitudinal expansion whereas decreasing transverse radii seem to 
indicate that also transverse expansion takes place (7]]^] . 

Resonances are known to be abundant in relativistic heavy ion collisions and they contribute significantly 
Qh! to pion production by, e.g., ^80% according to the Fritiof modebj^] Resonances affect the source such that 
■ it seems to have a larger life time or larger outward source size El]. The very long lived resonances are not 
resolved and reduce correlations. However, the reduction is also expected if the source is partially coherent. 
^ . The purpose of this paper is to do a combined analysis of resonances and flow both longitudinal and 
k^J ' transverse, and to calculate the nontrivial interplay. For that purpose the the correlation function for an 
incoherent source with plane wave propagation is extended to include resonances. By expanding in small 
relative momenta between the identical particles, one can extract the source radii analytically for general 
cylindrical symmetric sources with longitudinal and transverse flow. The effect of resonances and flow on 
the correlation function will be discussed and compared to recent NA44 data, which can extract sidewards, 
outwards and longitudinal source sizes. 
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II. SOURCE SIZE FLUCTUATIONS 



Assuming an incoherent source and plane wave propagation, the correlation function reduces to the simple 
Fourier transform ||r^ , |l3f 

C(q,K)_l± {Jd4xS{X)KW , (1) 

for two identical bosons/fcrmions (±) of relative and total momenta q — k\~ki and K = k\ + hi respectively, 
when q <C K. The source distribution S(x, K) — dN / cfixdVL describes the particle production at space-time 
point x = (t,x,y, z) of a particle with momentum k\ ~ — K. 

Resonances will here be included by assuming classical propagation. The additional distance travelled by 
the resonance, Ax = u^-Ar, where At is the life time of the resonance and u r its velocity, gives an extra 
phase, Ax ■ q, in Eq. (yj). When the life times are exponentially distributed with decay time r r = 1/T r , the 
correlation function becomes (see also El]) 
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Now, x refers to the space-time production point of the resonance. 

It is convenient notation to introduce the space-time integration and summation over resonances of an 
operator O by the average 

- = \Z r Jd*xS r (x,K)d\ 

{U) - \j:jd^s r (x,K)\ ■ (6) 

Likewise the fluctuation of that operator is 

a{d) = {d 2 )-{6) 2 . (4) 

The correlation function for bosons is then 

C(q,K) = l±{e iq - x {l-iq-u r T r )- 1 ) 2 . (5) 

Whereas the full g-dependence of the correlation function is generally quite cumbersome to calculate, the small 
q-dependence can be calculated directly by expanding (|^) to second order in q 

C(q, K) = 1 ± [l — a(q ■ [x + u r T r ]) - {[q ■ u r T r f) - 0((q ■ x) 4 )] . (6) 

The expansion is valid when R^q^ <d 1, where are the space-time extension of the system. If the source 
and therefore also the correlation function is approximately gaussian, then the leading term is sufficient to 
determine its shape. Notice that the resonance life-times contributes to the correlation function in Eq. (||) 
through the fluctuations in space and time, a{q ■ [x + u r r r ]), as well as contributing by the explicit term 
([q ■ u r T r } 2 }. 

Experimentalist often parametrize their pion and kaon correlation function by 

C(q, K) = 1 + A exp[-q*R 2 - q 2 R 2 - q 2 R 2 - 2q oqi R 2 ol ] , (7) 

where q , s ,z are the usual outwards, sidewards and longitudinal projections respectively of the relative momen- 
tum in a cartesian coordinate system where q; lies along the beam (or longitudinal) axis and q s is perpendicular 
to both the beam axis and K. The cross term with coefficient R i has recently been advocated by Heinz et al. 
H . Indications of this kind of asymmetry have been seen in experiments J6| . The prefactor A is found to be 
approximately unity for kaons and about half for pions in most experiments. The reduced pion correlations 
may be due to "coherence" in the source or to long lived resonances as will be discussed below. 

Because the momentum of the pair K breaks cylindrical symmetry and, with the beam line, determines the 
x-z plane, the average values for the z and x coordinates are generally non-vanishing in a longitudinally and 
transversally expanding system. Reflection symmetry in the y coordinate leads to (y) = and we can assume 
that the resonance velocity u r = (ut, u x , u y , u z ) has vanishing y-component, (u y ) — 0. Consequently, all cross 
terms with R s vanish as in (Q). Since q = (/3 q + Piqi,q , q s , qi), where (3i — Ki/Ka, we obtain from (Q) and 

(0) 

Rl = *{y) (8) 

R 2 Q = ct(x + u x t t - (3 (t + u t T r )) + {[(u x - l3 u t )T r ] 2 ) (9) 
R 2 = a(z + U z T r - Pl(t + U t T r )) + ([(u z - filUt)T r ] 2 ) (10) 

R 2 ol = (i X + U xTr - Po{t + U t T r )][z + U z T r - (3l(t + U t T r )]) 

- (x + u x r r - (3 {t + u t T r ))(z + u z r r - pi{t + u t T r )) + {(u x - (3 u t ){u z - Piu t )r 2 ) . (11) 

To evaluate the source radii, we need to specify the source further. Since the particle spectra have approxi- 
mately thermal transverse momentum distributions in the energy and p± regions considered here, we will here 
assume that all resonances are thermally distributed locally 

S r (x, K) = f r exp[-K ■ u{x)/T r ] p{x) . (12) 

Here, f r determines the fraction of particles coming from resonance r, T r is the "temperature" or more 
accurately the inverse of the m± slope of particles produced through that resonance, u(x) is the flow velocity 
at a given point in space and time, and p is the spatial and temporal source function. For longitudinally 
expanding systems the proper time r = \Jt 2 — z 2 and °q—\ ln((i + z)/(t — z)) are convenient variables whereby 
d 4 x = Tdrdrjclxdy. 

The source function is commonly parametrized by factorizing gaussians common for all resonances 
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p(T,T),x,y) exp 

T 

However, the detailed form will not be needed in the present analysis. Only the mean and the fluctuations in 
the quantities, e.g., (x) and c(x) = (x 2 ) — (x) 2 , are needed to leading order. When the experimental accurracy 
improves one may be able to measure finer details of correlation functions which then can determine higher 
moments. Whereas we can always translate our coordinate system so that (x) = (y) = (77) = 0, the average 
emission or freeze-out time (r) does not vanish and it determines the longitudinal extension of the source. 
The source in Eq. (Jl^) assumes cylindrical symmetry, i.e. c(x) = c(y). If one is able to determine the plane 
of reaction, that symmetry is broken and one may determine er(x) and cr(y) separately. 
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III. LONGITUDINAL EXPANDING SYSTEM 



Let us first analyse the situation with longitudinal expansion as in the Bjorken model, 
(cosh 77, 0, 0, sinh7j). Defining the mean rapidity as Y = (yi + 2/2V2, we find 



i.e. 



K ■ 



7tij_ cosh(?7 — Y) . 



(14) 



The resonance velocities are approximated by their local average over all directions which is simply the local 
flow velocity, i.e., u r — u. It is now straight forward to evaluate the source radii of d^[TT|). However, the 
gaussian in rapidity has the effect of moving the average rapidity, (77) (referred to as the "saddle point" in 
Refs. |§fr|), away from the mean rapidity, Y, of the two interfering particles. In the limit when a n 3> T/m± 
so that (77) ~ Y, the results simplify to 



^ 2 = *(y) 

R 2 = cr(x) + p x [<t{t) + <j(r r ) + (r 2 ) + (r + r r ) 2 tanh 2 Ya(r))] 
Rf = <t(t/)(t 2 + 2tt t + 2r 2 ) cosfT 2 (Y") 



-0±e(ji)(T + T rY sinh(F) cosli (F) 



(15) 
(16) 
(17) 
(18) 



where f3j_ = pj_/fn±. We observe that the source sizes are given in terms of the dispersion of the source, 
namely the fluctuations in spatial coordinates (<r(x) and c(y)), temporally (<t(t)), rapidity (a(j))), and in 
resonance life times (<7(r r )). For the simplified source of Eq. ( |13| ) the spatial fluctuations transversally are 
simply (j(x) = <r(y) = <r±, and the temporal fluctuations are a(r) = <J T . 

The fluctuation in rapidity consist not only of the spatial rapidity fluctuation in the distribution ( |i"3| ) but 
is strongly reduced by the additional rapidity dependence present in the thermal factor of (112) 



(7(77) 



/ ■mi 

V rp 



(19) 



However, as discussed above, it was assumed that a v T/mj_ in deriving Eqs. (p^-^8|) and therefore the first 
term in (|h]) dominates. When the various resonances have different temperatures, the average over T r will 
produce fluctuations in temperature, a(T r ), besides the average value (T r ). 

TABLE I. Resonances, that decay into pions and contribute to the pion correlation function. The fractions are pre- 
dicted in the Fritiof || and RQMD j[o| models for central S+Pb collisions at energies around 200 GeV/A, midrapidity 
and all mi, 



Resonance 


Decay 


m r 


Width 


CT r 


j Fritiof 


fRQMD 
Jr 


direct n 




140 MeV 




fm 


0.19 


0.33 


P 


7T7T 


770 MeV 


153 MeV 


1.3 fm 


0.40 


0.26 


A 


Nty 


1232 MeV 


115 MeV 


1.7 fm 


0.06 


0.12 


K* 


Kit 


892 MeV 


50 MeV 


3.9 fm 


0.09 


0.07 


E* 


Ett 


1385 MeV 


36 MeV 


5.5 fm 


0.01 


0.02 


LJ 


7r7T7r 


783 MeV 


8.4 MeV 


23.4 fm 


0.16 


0.07 


t 

n 


rjmr 


958 MeV 


0.24MeV 


821 fm 


0.02 


0.02 


V 


7r7T7r 


549 MeV 


1.1 keV 


1.2 A 


0.04 


0.03 


K° s 


7T7T 


498 MeV 


~ 


2.7 cm 


0.03 


0.07 


E,E 


n7r, nn 


1193 MeV 


~ 


4.4 cm 


0.00 


0.01 
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The resonances contribute to fluctuations by cr(r r ). When the temperatures or m^-slopes are the same for 
all resonances, the fluctuation in the resonance life times is 

<T r ) = <T r 2 ) - (r r ) 2 = frT 2 r ~ E frTr) 2 , (20) 
r r 

where f r is the fraction of pions arising from resonance r as, for example, given in Table I for the Fritiof 
|H and RQMD |hJ models. The resonance contributions may vary significantly with rapidity and transverse 
mass, i.e. the T r 's are different, and therefore the fraction at the relevant rapidity and transverse mass of the 
pair, f r (Y,m±), should be inserted in Eq. (pp|). 

In relativistic heavy ion collisions numerous resonances contribute to pion production as illustrated in Table 
I. The relative contributions from the various resonances have unfortunately not been measured very accurately 
and different models give a variety of results. This is illustrated by the Fritiof and RQMD models in Table I, 
where the feed-down to ir + around mid-rapidity are given for central S + Pb collision at 200 GeV/A corrected 
for detection efficiency. Note that the rescatterings in RQMD lead to a strangeness enhancement, which at 
mid-rapidities is mainly felt as an enhancement of the long-lived Kg and to less degree of hyperons. Feed-down 
from 3*, A, (ft, ... contribute by a few per mille each and have not been included in Table I. The feed-down to 
kaons is mainly from the K* by ~50% in Fritiof and ~5% in a thermal model |l5|). 

Since the typical minimum relative momentum, that can be measured between pions in relativistic heavy ion 
collisions, is of order q m i n ~ 5 — lOMeV, the resonance phase is large, q ■ t t ^> 1, for the long lived resonances 
r = 77, r/', Kg, X, E, and therefore the small q employed in Eq. (g) expansion is not allowed. Instead the 
resonance factor vanishes, |1 — iq ■ u r T r | _1 <C 1, for these long lived resonances and their contribution to the 
correlation function can be neglected. This can be taken into account by simply excluding these long lived 
resonances in the sum over resonances and the correlation function in ( |f3| ) is reduced by a factor 

Kea = (1 — f V — frf ~ f K% - /s,S,..) 2 ■ (21) 

The u> resonance has a life-time such that qm.in.Tui ~ 1 and should therefore be treated as an intermediate case. 
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S+Pb 200 GeV/A (NA44) 



long lived resonances 




40 60 
q [MeV/c] 



FIG. 1. Resonance contributions to the nix correlation function as function of outwards relative momentum, q 

o 



text). Curves include successively direct pions, p, A, K* , E*, ui and 77 + rj + Kg + E + E + 
Experimental data for central S+Pb collisions at 200 GeV/A JE| is shown for comparison. 



see 

from RQMD (Table I). 
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To illustrate the effects of resonances, the two-pion correlation function is shown in Fig. 1 as function of 
the out-ward relative momentum. The correlation functions are calculated from Eq. (^|) with the gaussian 
source of Eq. (13), which Fourier transforms to \p\ 2 = exp[— q%cr±]. The transverse radius y/a±_ = 4fm is 



assumed and the resonance fractions and life times arc from RQMD as given in Table I. The curves shows the 
contribution from resonances by adding them successively in the order of increasing life-times. It is assumed 
that the source is completely incoherent and that Coulomb effects have been corrected for. The long lived 
resonances lead to a smaller A value and the resonances in particular the u> makes the resulting correlation 
function steeper than a simple gaussian at small q a . Many experiments actually find that the pion correlation 
data is fit better by steeper functions (e.g., exponentials) than gaussians. 

For comparison, the ir + ir + correlation function measured in S+Pb collisions at 200 GeV/A [5j is also shown 
in Fig. 1. The data are collected for q s , qi < q cut = 20MeV/c ||. The cuts lead to a further reduction 

X = C(q ,q s < q cu t, qi < q cu t) — 1 = X res X cut , (22) 

where X cut = (1 — i? 2 g 2 ut /3)(l ~ Rfq 2 ut /3) when q cu tR s ,i^, and the quadratic expansion of Eq. (||) is valid. 
For the NA44 source sizes R s = 4.2fm and R t = 4.7fm and q cut = 20MeV/c one finds X cut = 0.87. In the 
comparison with experimental data of Fig. 1, the theoretical curves include the 13% reduction from X cu t- 



The data prefers more long lived resonances than predicted by Fritiof (see also |16|) and RQMD and there is 



no indication of the u> resonances. Due to strangeness enhancement of in particular Kg, the RQMD model 
describes the data better than Fritiof. In the recent Pb+Pb experiments at CERN, the small ~ 0.3 
measured [T^] may be explained by the larger q cut = 30MeV and larger source sizes R , s ,i — 6 — 7fm, which 
lead to a smaller \ cu t- No X cut applies when invariant relative momenta, qi nv — \fq*~, are employed which 
explains why Xi nv — C(qi nv = 0) — 1 generally are found to be larger. 

The measured values for the sidewards and outwards radii are very similar in high energy nuclear collisions. 
For example, in S+Pb collisions @ R a = (4.02 ± 0.14)fm and R s = (4.15 ± 0.27)fm so that R 2 S - R 2 S = 
(-1.1 ± 3.4)fm 2 . In Eq. @ the resonances makes a difference of a(r r ) = 1.2fm 2 and (t 2 ) — 3.0fm 2 in 
the Fritiof model even without including w's or the long lived resonances. The transverse velocities for the 
measured pions are /3± ~ 0.6 — 0.9. Thus the short lived resonances lead to a significant difference between 
R 2 — R 2 that is already between one and two standard deviations larger than the experimental value. This 
leaves very little room for additional fluctuations in source emission time, <t(t), and therefore the pions must 
be produced in a "flash" . 

Due to longitudinal expansion, both the expansion time of the source (r) and the lifetime of the resonances 
r r as well as their fluctuations contribute to R , Ri and R i. The longitudinal source size scales with rapidity 
and transverse mass as Ri oc (i/m_i_ cosh(F)) -1 as found in [Q and according to Eq. ( |l7j ) when er r) <C m±/T. 
Experiments confirm both the rapidity dependence ||[l7) and the m± dependence (see 0] and Fig. 2). If, 
for example, the thermal factor in ( p"2| ) is replaced by that of free streaming (~ exp(— mj_/T), see, e.g., Refs. 
^9|), then the resulting rapidity fluctuations would only be the second term in (|l9|). This does not produce 
the observed y and m± dependence which indicates that some thermalization must have taken place in high 
energy nuclear collisions before freeze-out. 

The source radii above can be Lorenzt boosted longitudinally to any frame by simply boosting the rapidity 
Y . In the longitudinal center-of-mass system, defined such that every pair is boosted to its c.m.s. system, 
the radii are obtained from the above by setting Y = 0. Thus the asymmetric term R D i vanish in Eq. ( |T§| ) to 
quadratic order but as pointed out in |8J terms of fourth order are non- vanishing. 

Most resonances produced in relativistic nuclear collisions cannot decay into two positively or negatively 
charged pions. Also, a single string cannot produce two particles with the same nonvanishing charge next to 
each other and an anticorrelation appears. These anti-correlations in a single resonance or string are, however, 
completely washed out by the abundance of resonances in relativistic nuclear collisions. 



IV. TRANSVERSE EXPANSION 



Besides longitudinal expansion the source may also expand transversally. Transverse flow has been studied 
with renewed interest |?],|8| since it decreases the source sizes R Q and R s with increasing mi, an effect that 
has recently been seen experimentally ||. The analysis in |7],|8| will here be extended to include resonances. 

A source expanding transversally with velocity v(x) has flow velocity 

u = j(v)(cosh(r)),v x (x),v y (x),$mh(r])) , (23) 

which gives 

K ■ u = m_L7(u)(cosh(?? - Y) - P±v x ) . (24) 
Assuming the source is expanding as v x = vxj \joi_ and similarly for v y one finds when d 2 «1 
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1 + 2„,>f + ^/ 
1 + v z m±/l r 

R 2 = R 2 S + (Pi + v 2 )(t 2 ) + Pl [(1 + v 2 )a(r r ) + <j(t) + (r + (1 + « 2 )r r ) 2 tanh 2 Ftrfa)] (26) 
= <K»?) [( t2 > + 2(1 + w 2 )(rr r ) + 2(1 + « 2 ) 2 (r 2 )] cosh^Y) (27) 
i?2 ; = -0 ±a (ji){T + (1 + v 2 )r r ) 2 sinh(F) cosbr 2 (F) . (28) 

One notices the factor (1 + v 2 m_\_/T)~ 1 which was found in earlier analyses |7],|| and which leads to smaller 
apparent sources at large m± . Resonances do not change this dependence but lead to generally larger apparent 
sources by additional terms of order vr r /a±. 

The transverse flow results in smaller transverse momentum slopes or equivalently larger apparent or effective 
temperatures, T e ff. For a fixed m± slope, the local temperatures are "red shifted" approximately as T — 
T e ffy/ (1 — v)/(l + v) |Q as the transverse flow increases. Consequently, the transverse flow results in smaller 
source sizes at large m± not only due to the explicit v 2 terms in Eqs. (|25|-^8|) but also by reducing the 
temperature. 
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200 300 400 500 600 

m t [MeV] 

FIG. 2. Effects of transverse flow on the sidewards radius as function of transverse mass (see text). Curves correspond 
from top and downwards to v/c = 0, 0.2, 0.4, 0.6, 0.8 . Experimental 200 GeV/A S+Pb data || on R s , R , and Ri are 
shown (with squares, circles and diamonds respectively) for pions at two different m± and for kaons at m± — 550MeV. 

In Fig. 2 the m± dependence of the sidewards radius is shown for T e ff — 170MeV ||, a± = 5fm, and for 
various transverse flow velocities. As seen from Eq. (|25|), the scale is mainly determined by aj_ whereas the 
curvature depends on the transverse flow velocity. The decreasing size with increasing m± can be reproduced 
with flow velocities of order v/c ~ 0.6 — 0.8. These transverse velocities are smaller than those obtained 
by Csorgo and Lorstad. Resonance fractions from RQMD (Table I) have been employed but using fractions 
from Fritiof does not change the dependence on m± by much. The temperature is assumed to decrease with 
transverse flow velocity according to the red-shift formula. Including resonance fractions varying with m± as, 
e.g., in RQMD makes only minor differences in the m± dependence of the source sizes. 
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DISCUSSION 



A general description of correlations and source size dependence on resonances and flow has been given. By 
expanding to second order in relative momenta the source radii were extracted analytically and a nontrivial 
interplay between the effect of resonances and flow was found. The source sizes receives contributions from 
fluctuations in transverse spatial directions, source life time, rapidity and resonances life times. The long 
lived resonances lead to a reduction in the correlation function and can explain most of the reduction in the 
correlations function measured experimentally, A 7 ™ ~ 0.5 — 0.6, when the experimental cuts are included. 
Due to strangeness enhancement in RQMD of in particular the long lived Kg, it gives a smaller A in .better 
agreement with experiment than Fritiof. However, the similar values for A 71 " 71 " in pp and pp (A ~ 0.40 §fel) and 
p-nucleus (X pA ~ 0.41 |J) and nuclear collisions do not indicate that rescattering, thermalization, strangeness 
enhancement or other collective effects in nuclear collisions are affecting the resonance production and A 
significantly. There is no indication of the w-resonance in the data, which is otherwise abundant in most 
models. From measurements of the KK correlation functions one finds \ KK = 0.83 ± 0.08 ||), which is close 
to A CU f, and therefore no long lived resonances decaying into kaons are required. This is predicted in the 
models discussed above where kaons mainly receive feed-down from the relatively short-lived K* . 

Even though resonances and cuts can explain most of the reduction in X n7T found in nuclear collisions, 
~ 0.5 — 0.6 |j 21 14|], it does not exclude a partially coherent pion source. However, the kaon sources 



measured, \ KK /\ cu t ~ 1, are apparently incoherent. Furthermore, if particles are produced coherently within 
a coherence length, £coh> m space and time, the reduction in the correlation function would be of order 
A ~ 1 — (^coh/R) 4 '. The similar A values found in pp, pp, p-nucleus and nucleus- nucleus collisions would then 
require that the correlation length scales approximately with the source size. 

Other effects, that might explain the low A 7r7r , are final state interactions and Coulomb repulsion. Strong 
interactions are generally found to be insignificant in nuclear collisions due to their short interaction range as 
compared to nuclear length scales. Coulomb repulsion reduce the correlation function significantly at small 
relative momenta, which is usually corrected for by the Gamov factor T = ij/(e 7, — l), where fj — 2irme 2 /Q inv ~ 
6MeV/Qi„„ for pions. Coulomb screening effects are only of order a few percent p2| but as they reduce the 
Coulomb interactions, they lead to a smaller A. 

Experimental data seem to indicate both longitudinal and transverse flow. Longitudinal flow as in the 
Bjorken model seems to agree with the rapidity dependence of the longitudinal source radius. The decreasing 
transverse source sizes with increasing transverse mass of the particles can be explained in a simple model 
with transverse flow of order v = 0.5 — 0.8c at the surface when resonances are included. This assumes an 
apparent blue-shifted temperature that is fixed by the transverse momentum slopes of particle spectra. The 
source sizes have been measured for two values for the pion transverse mass and one kaon transverse mass. 
The smaller kaon radii can thus be explained solely by transverse flow. The conventional explanation is that 
the smaller kaon scattering cross section results in earlier freeze-out at smaller radii than for pions. However, 
such a cascade picture leads to an extended particle emission time typically of the same order as the expansion 
time, <t(t) ~ (r), and not a sharp freeze-out in time as is indicated by the similarity of the sidewards and 
outwards radii. Also the freeze-out in hydrodynamic models at a constant temperature takes about as long 
as the expansion time. This is in contradiction with the very similar outwards and sidewards transverse radii 
measured in particular when the resonance life-times are taken into account. 

Clearly, more data is crucial in order to describe the sources and determine their size, lifetime, longitudinal 
and transverse expansion, so that we can discriminate between these models. In particular, it would be most 
useful if one could measure the feed-down from resonances independently. One would then be able to predict 
the non-gaussian correlation function (see Fig. 1) and separate the effect of resonances from coherence, final 
state interactions etc. Only when resonances are properly accounted for and their effect on the correlation 
function, A, the source radii and life time are understood and corrected for, can we extract the "bare" source 
sizes and the flow reminiscent of the earlier phase of hot and dense matter created in high energy nuclear 
collisions. 
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x The pp and pp data are, however, at higher collision energy, %/s = 63GeV, and at larger relative momentum, 
q > 50MeV/c. 
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